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ABSTRACT 

Ensemble sets of simulation experiments were conducted with a single column model (SCM) 
using the Goddard GEOS II GCM physics containing a recent version of the Cumulus Scheme 
(McRAS) and a biosphere based land-fluxes scheme (SSiB). The study used the 18 July to 5 
August 1995 ARM-CART (Atmospheric Radiation Measurement-Cloud Atmospheric Radiation 
Test-bed) data, which was collected at the ARM-CART site in the mid-western United States and 
analyzed for single column modeling (SCM) studies. The new findings affirm the earlier 
findings that the vegetation, which increases the solar energy absorption at the surface together 
with soil and soil-moisture dependent processes, which modulate the surface, fluxes (particularly 
evapotranspiration) together help to increase the local rainfall. In addition, the results also show 
that for the particular study period roughly 50% of the increased evaporation over the ARM- 
CART site would be converted into rainfall with the Column, while the remainder would be 
advectedout to the large-scale. Notwithstanding the limitations of only one-way interaction (i.e., 
the large-scale influencing the regional physics and not vice versa), the current SCM simulations 
show a very robust relationship. The evaporation-precipitation relationship turns out to be 
independent of the soil types, and soil moisture; however, it is weakly dependent on the 
vegetation cover because of its surface-albedo effect. Clearly, these inferences are prone to 
weaknesses of the SCM physics, the assumptions of the large-scale being unaffected by grid- 

scale (SCM-scale) changes in moist processes, and other limitations of the evaluation 
procedures. 
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1. Introduction 


Earth’s biosphere can influence fluxes in land-atmosphere interaction through four primary 
controls, namely by: (i) increased absorption of solar radiation within leaf organizations of 
vegetation canopies; (ii) increased evapotranspiration with access to root-zone soil moisture; 
(iii) stomatal control that generally stifles evapotranspiration during warm and/or dry 
episodes; and (iv) contributing surface roughness on the scale of turbulent eddies, which in 
turn increase the boundary layer depth and cross-isobaric moisture convergence in a typical 
surface-low pressure episode. In dry regions with small evapotranspiration, an increase in 
the surface-albedo reduces the solar absorption, land-surface fluxes, and rainfall (e.g., Sud 
and Molod, 1988; Dirmeyer and Shukla, 1994). Likewise, reduction of evapotranspiration 
and/or land-surface roughness -- a potential effect of removal of vegetation - also decreases 
rainfall (e.g., Sud and Smith, 1985). In fact, a discernible dependence of rainfall on 
vegetation has been simulated and/or discussed in several papers (e.g., Dickinson, 1 980; 
Anthes, 1984; Avissar, 1992; Sud and Fennessy, 1982, 1984). Moreover, the dependence of 
summer season precipitation on biospheric processes can be inferred from the first 
principles of moist-convection as enunciated by Arakawa and Schubert (1974). The very 

same principles have been used in the design of a whole generation of physically-based 
cumulus schemes, including FIAS (Moorthi and Suarez, 1 992), and McFtAS (Sud and Walker, 
1999). Sud et al. (1 993, 1995) pedagogically argued and demonstrated that each vegetation 

process, when limited to influence only the column-atmosphere, helps to promote moist 
convection. Specifically, a surface albedo increase leads to near-surface energy deficit, 
reduces CAPE (Convective Available Potential Energy), and suppresses moist convection. Ch 
the other hand, simulations with altered proportions of evapotranspiration and sensible 
fluxes, for scenarios in which the sum of the surface fluxes remains invariant, reveal that 
evapotranspiration is a more desirable surface-flux for promoting convective rainfall 
because it enables CAPE-accumulation before "turning on" the moist-convection. Larger 
CAPE naturally enables convection to reach deeper into the atmosphere thereby producing 
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more precipitation and associated atmospheric warming. Since convective (as opposed to 
stratiform precipitation) produces less fractional cloud-cover, it naturally allows more 
insolation to reach the surface of the Earth, which further boosts the energy needed for 
surface fluxes. Together with the observational analyses (e.g., Otterman et al . , 1 990; 
Nicholson, 1985; and Skole and Tucker, 1993), such findings provide a rational physical 
basis for understanding the influence of vegetation on precipitation. 

Even in several earlier studies, a warmer and drier climate was simulated in response to 
Amazonian deforestation (e.g., Dickinson and Henderson-Sellers, 1988; and Henderson- 
Sellers et al., 1993), even though there were significant differences in the simulated 
rainfall. The outcome became more perplexing when results of several GCMs (each with its 
own interactive biosphere) were examined side-by-side. Table 1 of Hahmann and Dickinson 
(1997) contains a comparison of evapotranspiration, moisture convergence, and local 
rainfall simulated by fifteen different GCMs used for simulating the influence of Amazonian 
deforestation. Even though a majority of models simulated a substantial decrease in rainfall i n 
response to deforestation, a few models simulated an increase. This was caused by the 
competing effects of moisture convergence produced by the thermal heating of the lower 
troposphere vis-a-vis the moisture deficit produced by reduced evapotranspiration in 
deforestation (Eltahir and Bras, 1993). Recognizing that some differences in the Hahmann 
and Dickinson (1997) comparison would be due to differences in the location and the extent of 
the imposed deforestation anomalies or in the simulated climatic index of dryness (Koster and 
Suarez, 1999), it is also abundantly clear that the simulated findings continues to be model- 
dependent. Since models naturally simulate the synthesized outcome of modelled interactions 
between physics and dynamics, such variations in results point to the lack of satisfactory 
representation of modelled processes and their interaction with dynamics. We submit that 
there is a fundamental need for: (i) more extensive validation of modelled processes in GCMs; 
and (ii) devising more ingenious ways of examining the land-atmosphere interaction problem 
to arrive at a better answer. Consequently, we begin by examining the evidence of our model’s 
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performance in a Single Column Model (SCM) using the prescribed large-scale forcing. Only 
after the evaluation are judged to be satisfactory, will we use the SCM to analyze its response 
to biosphere-precipitation interactions? 

We examine the influence of the biosphere on the summer precipitation over the ARM-CART 
site of the mid-western Great Plains (for site details, see Randall and Cripe, 1999 Fig. 1 , 
page 24,531). This region is affected by strong moisture advection by the low-level jet, 
which is a manifestation of the large-scale forcing. Our Single Column Model (SCM) of the 
atmosphere is coupled with simple SiB (SSiB) for simulating biosphere-atmosphere 
interactions that in turn are coupled to a highly resolved 1-D soil hydrology model that will 
be shown to yield a very accurate solution of Richard’s equation. A state-of-the-art cloud 
scheme, called McRAS (Sud and Walker, 1999), is already coupled to the GEOS II GCM physics 
and has been shown to produce realistic results. The goal is to simulate most realistic soil- 
vegetation-atmosphere interactions, with the lower boundary condition at 5m below (water 
table may be higher) the surface of the soil. Assuming that such a model is realistic for the 
Great Plains region, what can we learn about the behavior of the column atmosphere to 
changes in surface fluxes? The question we pose is: if the large-scale atmosphere was 
essentially unaffected by much smaller-scale changes in the land-surface fluxes due to 
changes in vegetation cover or soil type or soil moisture, how would the vegetation-soil- 
rainfall interaction be affected locally? To examine this, we force the SCM to interact with a 
suite of soil types, soil moistures, and/or vegetation covers fractions. The July 1995 ARM- 
CART data represents wet conditions and we only perturb it towards drier conditions. We will 
describe our SCM, the cloud scheme called McRAS, and the new 100-layer soil-hydrology 
model in section 2. We also show concrete evidence of the intrinsic reliability of the SCM’s 
physics and make a case for the suitability of this model for the proposed investigation in 
section 3. We then show key results of SCM studies with each of the five soil types adopted i n 
the ISLSCP Initiative I data with a full range of soil moisture and vegetation covers in section 
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4. Finally we will discuss the implication of the findings for ecosystem and biosphere 
managers in section 5. 
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2. Single Column Model 


Single column models (SCMs) have a long history of use for evaluating numerical models 
representing boundary-layer processes (e.g., Hoffert and Sud, 1976), land surface hydrology 
(Koster and Eagleson, 1990) and moist processes including convection and downdrafts (e.g., 
Cheng, 1989). Randall et al. (1996) discussed a variety of innovative uses of SCMs and 
illustrated how a single column model (SCM) would be a useful device for model 
parameterization evaluation and/or scientific hypothesis validation. We will use our SCM to 
validate the hypothesis that land-surface biospheric processes, particularly 
evapotranspiration, have a positive feedback effect on convective precipitation. To meet these 
goals, we force the SCM with observational data. This enables the investigator to differentiate 
between the mathematical and/or numerical characteristics of a particular parameterization 
vis-a-vis the model-simulated reality of the atmosphere. Our earlier results were that all 
biospheric processes help to promote the local rainfall, particularly the convective type 
everywhere, and that the inference is expected to hold for the mid-western United States 
typified by the ARM-CART site. However, the magnitude of this feedback was not discerned i n 
those studies (sud et al., 1993). 

The particular SCM is extensively described in Sud and Walker (1993) who used it for 
evaluating convective downdrafts using GATE Phase III data. More recently, Sud and Walker 
(1999) used the same SCM for evaluating the performance of a new Cloud microphysics with 
Relaxed Arakawa-Schubert Scheme (McRAS). The ARM-CART data is a new dataset for SCM 
use. Recently, it was analyzed and used for evaluating as well as model intercomparison of 
mid-latitude cloud processes (Ghan et al., 2000). Randall and Cripe (1999) discuss three 
different ways in which an SCM can be used for testing physical parameterizations. In the 
very basic form, i.e., without flux adjustments or artificial relaxation, which the authors call 
“Revealed Forcing", the specific humidity, q, tendency of the column atmosphere can be 
expressed by: 
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Here V • Vq is the domain-averaged horizontal, and — ^ - is the domain-averaged vertical, 

3p 

flux divergence tendencies of q. Clearly, these tendencies are in the flux form and contain the 
sum of advective plus divergence tendencies. All the modelled physical processes are contained 
in the P-tendency, the last term of RHS. Equation (1) can be recast with the specific 
humidity q split three ways- (i) into observed, q obs , (ii) land-surface fluxes forced anomaly, 
q veg , and (iii) SCM- errors forced effect, q me , as follows: 
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The observed tendency can now is transformed into the basic advective form, which merely 
alters the first term of RHS to give: 
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A simulation, with no vegetation or soil-type anomalies, will have only the first and third 
terms of (3b) contained in square brackets of RHS, plus its tendency due to physical 
adjustment, P conl . The land-surface fluxes anomaly simulation will have all the three terms, 
but each with its own physical adjustment term, P veg . By appropriately accounting for these 
terms, we can delineate the vegetation and/or soil anomaly influence on the ensuing 
precipitation. 
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First let us rewrite (3b) in the following form: 
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where [O A T] is the observed advective tendency. The control and the anomalous land-surface 
fluxes simulations can be performed with the SCM using observed advective tendencies from 
the ARM-CART observations that get modulated by the second term of the RHS, which 
comprises of model errors and vegetation plus soil-type anomaly influences. The model is 
assumed to be realistic, but the model errors can be reduced largely, if not entirely, very 
much the way GCM-errors are expected to do in producing anomaly minus control differences. 
This would happen precisely if the errors had a linear influence. By subtracting simulation 
(4a) from (4b) or by comparing them side-by-side, we can hope to infer vegetation- 
produced tendencies with resultant physical interaction effects that would appear in a GCM 
with full hydrodynamics. Our analysis technique is a modification of the artificial relaxation 
methodology of Randall and Cripe (1999), because the advective equations have been recast to 
delineate the influence of surface fluxes on the ensuing changes in the column-atmosphere. 

The second term in the square brackets of equations (4 a & b) represents the horizontal and 
vertical advective tendencies that are specific to the control and vegetation and/or soil 
anomaly simulations. To estimate them, we use the grid-cell domain-averaged horizontal 
wind, V, and vertical pressure velocity, to, from observations. Both of them are provided in 
the ARM-CART SCM data. Our mode of testing would be the “revealed forcing” mode of Randall 

and Cripe (1999) with aforestated corrections. The second term of (4a) contains the 
following: 

^9 h or _ Yobs *^(qmod ~ qobs) / 

a t ~ Xx " (q m od -qobs) v *v obs , (5) 
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where AX is the horizontal-length scale of the domain, while V obs is taken to be constant. Here 
qmod generically refers to the modelled q, which is different from q obs due to model errors 
and/or vegetation anomaly influences, as stated before. 

The second part is for the vertical flux divergence of q, or potential temperature, 0, as the 
case may be. For the moisture, its form is: 

dq vert _ _ ^ ^obs (^mod ~ 9obs) , _ . 

at a p ( 6 } 

Here <o obs is a function of p and must be multiplied with q-anomaly to conserve the anomalous 
moisture flux in the entire column during vertical advection. If we run all the cases involving 
different vegetation covers and soil-types by invoking (5 & 6) with or without including 
vegetation/soil effects, the influence of vegetation and/or soil-type on the parameterized 
physical interactions and observed advective tendency of q will be revealed. Whenever q of the 
SCM is different from the observed, the q advected in the horizontal and vertical direction is 
adjusted. The formulae are conceptually similar to that of Randall and Cripe (1999), except 
that our horizontal time constant T is simply: 

T hor = Horizontal Grid-length scale/Horizontal wind velocity. ( 7 ) 

A parallel equation can be written for other tendencies of which potential temperature ©is 
directly relevant. If the SCM’s q or © are influenced by surface processes, there is an 
appropriate adjustment for all material quantities and conserved variables. 

2.1 Cloud Model 

The cloud model of our SCM is McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert 
Scheme) and was designed and developed with the aim of simulating realistic moist processes, 
clouds, and cloud-radiation interactions in GCMs. McRAS distinguishes and parameterizes three 
types of clouds: convective, stratiform, and boundary layer as shown schematically in Figure 1 . 
The Relaxed Arakawa-Schubert scheme (RAS) of Moorthi and Suarez (1992) generates 



convective clouds, but the cloud-microphysics follows Sud and Walker (1999). The simulated 
convective clouds transform and merge into stratiform clouds on 1-h time-scale, while the 
boundary-layer clouds merge into the stratiform clouds instantly. The cloud condensate converts 
into precipitation following the auto-conversion equations of Sundqvist (1988, 1993) that contain 
a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of 
cloud condensate by precipitation. All clouds convect, advect and diffuse vertically (but not 
horizontally in the SCM-mode of testing) with a fully interactive cloud-microphysics. In this 
way, the life cycle of the cloud is related to cloud dynamics and microphysics, while their optical 
properties are derived from the statistical distribution of hydrometeors and idealized geometry of 
clouds. An evaluation of McRAS in a Single Column Model (SCM) with the GATE Phase III data 
was given in Sud and Walker (1999). It showed that together with the rest of the model 
physics, McRAS simulates the observed temperature, humidity, and precipitation without 
discernible systematic errors. The time-history and time-mean incloud water and ice 
distribution, fractional cloudiness, cloud optical thickness, origin of precipitation in the 
convective anvils and towers, and the convective updraft and downdraft velocities and mass 
fluxes all simulate a realistic behavior, even though all the diagnostics were not verifiable 
from the available data. 

In order to determine the suitability of McRAS for ARM-CART data use, McRAS was evaluated 
for with two ARM-CART datasets: CASE 1 and 3. The cloud scheme was run using data without 
the model error adjustments discussed above. The results are shown in Figs. 2a through 2c. 
Even though the model with GATE PHASE III data (Sud and Walker, 1999) was virtually 
error-free, we noted a few systematic errors in this evaluation. The SCM produced somewhat 
less rainfall than observed, and its lapse rate near the surface was very different from the one 
given in the data. For now, we must contend with the former because there are negligible 
column integrated humidity errors in these simulations, while the horizontal advective 
tendency, and the surface evaporation data were prescribed from observations. Therefore, 
any significant mismatch between the observed and SCM-simulated time mean rainfall reveals 



a lack of moisture conservation in the forcing ARM-CART data because our SCM naturally 
employs an internally consistent moisture conserving schemes. 

2.2 Land Model 

Our basic land model is the Simple SiB (SSiB) due to Xue et al (1991). This model has been 
tested with several datasets, but due to factors discussed in the Section I, land models of the 
present-day need to be constantly evaluated and improved, even region by region. A brief 
discussion of these improvements can be found in Mocko and Sud (2000, in preparation). Cne 
also notices SSiB parameterizes inter-layer hydraulic conduction among its layers with 
Richard's Equation with several arbitrary assumptions, which, in a crude way, help to give 
reasonable vertical fluxes. We felt there was a need to generate more accurate solutions of 
Richard’s equation with better numerical rigor. Consequently, we solved the equation with 
much higher resolution. We finally settled on 100-layers model with 5cm resolution (Fig 
3a). It is not suitable for GCM work but is useful for developing a parameterization of 
hydraulic conduction for land, which has no subgrid-scale variability. Indeed, lack of subgrid 
scale variability is an obvious limitation of the current SCM, but that is a separate issue and 
we reserve its treatment for a later time. Thus, the soil model has the bottom at 5m depth 
while the vegetation (root zone) layers of SSiB may vary regionally as a function of soil type 
and vegetation covers. SSiB was evaluated with ISLSCP Initiative I data under GSWP 
(Dirmeyer et al, 1999) as well as against other simple Land-Schemes (Mocko and Sud, 
1998). The dataset covers a period of twenty-four months: January 1, 1987 through 
December 31, 1988. In this data, only five (Table 1, Fig. 3b) soil types were used to 
represent the entire world, while there could be 10 different biomes; nevertheless, any one 
grid cell could have only one soil type and one biome. In other words, tiling is not feasible in 
the SSiB design. The subgrid-scale variability is addressed by allowing full flexibility in 
choice of fractional vegetation cover and leaf area indices. The new 100-layer model was 
extensively tested with ISLSCP Initiative I data without use of any tuning factors for a range of 
conditions-from extremely dry to highly moist. The model is able to generate realistic soil- 



water transport time-scales for all soil types and plausible vegetation processes. 

As an example, we show in Fig. 4a how altering the vegetation in the ARM-CART site from a 
desert to a fully vegetated ground cover would affect the soil moisture and land-hydrology. 
Desert conditions lead to less evaporation and higher soil moisture and soil water drainage or 
runoff, while full vegetation cover produces more evaporation and less soil moisture and 
runoff. These differences are clearly noted in the difference maps (Fig 4b). Such evaluations 
were made for monsoonal India, Sahel, Mississippi basin, and over a dozen other regions of the 
world (not shown). All these regions simulated very realistic behaviors. Without the aid of 
vegetation to remove the soil-water by transpiration, the soil moisture fraction in the month 
of June at 1-1.5 m depth can become as high as 0.6 versus 0.3 for the simulation with 
natural vegetation cover (top Panel, Fig. 4b). With 99% vegetation cover, there is a 
systematic reduction of soil moisture, particularly in the spring season in which there is 
considerable soil moisture uptake by transpiration. However, in late summer, the soil 

moisture does not go much below its value for the natural ground cover because vegetation 
starts to wilt in the dry periods. Similar results are noted for Sahel, which has an extended 
dry period (Fig. 5a & b). In fact, with the onset of the rainy season, the soil wetting time for 
the deep soil is of the order of 10-20 days as can be noted from the slope of the soil moisture 
isohalines going into the deeper regions of the soil. Once wetted, the adjustment time scale is 
much smaller. Regardless, wetter (drier) soil up to 3m depth is noted for minimal 
(maximal) vegetation cover simulations. Having satisfied ourselves of a reasonably realistic 
performance of the hydraulic model, we linked the 100-layer soil-SSiB model to the 
atmospheric SCM to conduct the vegetation-soil-atmospheric interaction studies. 


3. Design of the Experiment. 

The SCM (Sud and Walker, 1993 & 1999) with GE06 II GCM physics and McRAS now has a 
100-layer soil model (Koster et al. , 2000), which hydrologically responds to the natural 



vegetation in SSiB. To evaluate the influence of vegetation on the rainfall in the mid-western 
region, we employ ARM-CART data for Case 1. This is a 3-week dataset: July 18, 1995 
through August 5, 1995, averaged over a 300x300 Km region in Kansas and Oklahoma (called 
ARM-CART site). Even though the ISLSCP Initiative I data indicates roughly 60% vegetation 
cover, we allowed the vegetation cover to vary from 1% to 99% in this investigation. This 
helped to delineate the influence of vegetation on the simulated rainfall. Our experiment is the 
SCM equivalent of a regional experiment in which the large-scale forcing is prescribed 
through Analysis of Observations, whereas the region is free to respond with changes in its 
surface fluxes, but with parameterizations attempting to represent the changes in meso-scale 
circulation. The influence of vegetation on surface fluxes and therefore the temperature and 
humidity will be communicated in the vertical through physical processes, namely 
turbulence, dry convection, and moisture-cloud-radiative forcing of the column atmosphere. 
Clearly, we also have the freedom to alter the soil type with a better-resolved calculation of 
soil-moisture dependent soil properties. Since ISLSCP Initiative I recommended 5 possible 
soil types for the world, and there is a possibility of large variations in soil types in a typical 
300X300-km grid, we tested the SCM with each soil type using the same observed forcings. 
The initial soil moisture profile were obtained by making a long run with ISLSCP Initiative I 
data, and because of the differences in the hydraulic conduction properties of different soil to 
episodic rainfall events, the soil water profiles generated with different soils were so 
different in the structural complexity that we had to run each soil type and vegetation cover 
for six months prior to the simulation to get a reasonable estimate of the soil-moisture at the 
initial time. We ran a suite of 30 runs - one each for six values of fractional vegetation 
covers with each of the five soil-types. These are called control runs. The initial soil 

moisture profiles for each of the thirty cases are given in Fig. 6. We further produced an 
ensemble of five sets of thirty simulations each, in which the initial soil moisture profile of 
the 60% vegetation cover for each soil type was varied to introduce additional soil dryness 
0%, 25%, 50%, 75% and 100% of soil moisture deficit = 1- f s , where f s is the soil 
moisture. These runs were used to analyze the relationships among soil-moisture, soil-type, 



vegetation cover, and rainfall. 


4. Results. 

The objective of the current study is to determine the influence of vegetation on the rainfall in 
the mid-western United States. This test bed is a 300x300 Km region in Kansas and 
Oklahoma (called ARM-CART site). The site provides surface fluxes, top of the atmosphere 
radiative inputs, and advective and divergent heat, moisture and momentum tendencies for the 
single column model (SCM) atmosphere at 16 levels. The data were analyzed in 4-D VAR and 
made available at 3h intervals. We used the data to perform a 30 ensemble Control and 5 sets 
of ensembles of thirty cases each to examine the influence of soil moisture anomalies as 
described in section 3. The thirty soil moisture initializations were based on a 6-month 
simulation produced with ISLSCP Initiative I data for 1987 (a wet year that fits well with the 
wet summer of 1995) for each of the five soil types and six vegetation covers. Since the idea 
is to produce somewhat realistic soil moisture at the initial time, such a method is 
satisfactory for addressing the problem at hand. This is a compromise because the forcing 
data for the prior six months of 1995 was unavailable. We then performed the time series 
analysis and intercompare the simulations for each soil type and vegetation cover. The five 
sets of 30 runs each were produced in which the initial soil moisture profile varied by soil 
type only regardless of the vegetation fraction. The soil-moisture values were chosen to be the 
one for 60% vegetation cover (Fig. 6) for the first set. The soil moisture deficit (saturation 
value (=1.0) minus the initial soil moisture at each of the 100-soil layers) was increased in 
steps of 25% for each of the subsequent four sets making it 100% more for the last set. This 
generated a total of 150 simulations. They were used to address more fully the influence of 
initial soil moisture on evapotranspiration and its influence on precipitation. 

4.1 The Control Simulation 

Figures 7a-c show the temperature, relative humidity, and precipitation errors (simulated minus 



observed) at 3h intervals for soil type 3 and 60% vegetation cover — a combination is nearest to 
the observed land-surface conditions. All other diagnostics akin to those discussed in Sud and 
Walker (1999) were produced and seemed realistic, but only three of those that could be verified 
against observations are shown. A fully interactive simulation with the same SCM, but with 
surface fluxes prescribed from observations without any other relaxation, were noted to be 
reasonably realistic (Ghan et al. 2000). The current simulation is fully interactive and has SSiB 
coupled to a 100-layer subsoil hydrology model (Koster et al., 2000). Naturally, it generates its 
own surface fluxes. However, the atmospheric simulation contains an additional adjustment for 
model errors as well as vegetation anomalies as discussed in Section 2. With these adjustments, 
it yields even more realistic time-mean atmospheric soundings, i.e., time series of the vertical 
profiles of temperature, Fig. 7a, and humidity, Fig. 7b. In the 3-week averages, (shown on the 
far right), the simulated time mean biases are at an acceptable level, while the instantaneous 
biases keep fluctuation about the mean. We note that there are two episodes of about 4°C cooling 
with a preponderance of 1-2°C warm (cool) patches at around 100 (400) hPa level, respectively. 
In the moisture field as well, there were no systematic biases. However, one could infer that 
near the surface the model atmosphere is a little too warm, which compensates for the cooler 
environments atop. However, recognizing that there is virtually no systematically growth in 
temperature error, we contend that the simulation is quite reasonable. Since the simulated 
precipitation is less than the observed (Fig. 7c), it implies either the simulated evaporation is too 
low or the horizontal moisture flux divergence is too high. As we will show, the 
evapotranspiration is quite realistic for the 60% vegetation cover, therefore the prescribed 
horizontal moisture convergence may be the culprit, and this could be caused by a model bias 
and/or the input data. Because these errors do not relate to any systematic errors in the 
accompanying fields, it seems to indicate that the model is performing well within the limits of 
observational errors. 

The evapotranspiration and sensible heat flux time-series, shown in Fig. 8a & b, contain 
observations (black) and three vegetation cover cases with 1%, 60% and 99% vegetation (dashed 



lines) for soil type 3. Some differences between the observed and simulated surface fluxes can 
not be explained by vegetation variations, but recognizing that the SCM has its own clouds, 
water vapor, long and short wave radiation parameterizations, the simulated fluxes can be 
considered realistic particularly because they depict good diurnal amplitudes and phases. The 
only purpose of showing them here is to convey to the reader that the SCM has the ability to 
simulate realistic surface fluxes and respond to rainfall episodes that suppress surface fluxes by 
cloud shielding of solar radiation. During the continued dry period: July 26, 1995 through July 
30, 1995, the simulated sensible flux (evapotranspiration) is less (more) than the observed. This 
suggests that the parameterized vegetation had access to more soil moisture even though the 
simulated rainfall is lower. This is of some concern because it points to deficiency of the land 
model or biases in initial soil moisture profile (initialized with 1987 data), but the overall 
temporal structure of evapotranspiration is quite reasonable for addressing the current problem. 
Even though our model continues to produce reasonable precipitation time series for both cases 
(Fig. 8c) the current results show that the model can be considered satisfactory despite the 3- 
week mean precipitation being slightly less than the observed. 

4.2 The Anomaly Simulation 

As mentioned earlier, we used the SCM with above attributes to run 30 cases consisting of all 
possible combinations of 5 soil types (Fig 3b) and six discrete vegetation covers: 1, 20 40, 60, 80, 
99%, respectively. The initial soil moisture profiles were determined from a 6-month soil moisture 
initialization run using 1987 ISLSCP Initiative I data. This compromise allows each soil-type and 
vegetation-cover combination to use an initial soil moisture profile appropriate for it, allowing us to 
address the problem of the influence of soil type and vegetation cover on precipitation using 
relatively more realistic initial soil moisture profile. A total of 150 simulations were subsequently 
used to establish the robustness of our findings to the initial soil moisture using arbitrary soil 
moisture anomalies as discussed above. 


a) Evapotranspiration: 



Rgure 9a shows that the simulated evapotranspiration is a strong function of vegetation 
cover and soil type for the 3-week period. It is highest for soil type 5 and lowest for soil 
type I. The relative patterns are robust regardless of the initial soil moisture-whether they 
are taken from 6-month initialization for each vegetation cover or from initial soil moisture 
for 60% vegetation cover case. The hydrologic characteristtcs of different soil types are 

distinguished mainly by soil water suction potential, ig, and hydraulic conductivity, k, as a 
functional of soil moisture fraction through: 


and 


where Vs and k s are the corresponding saturated soil values; f w is the soil moisture fraction 
and B is Clapp and Homberger (1978) parameter. An increase in evapotranspiration for a 
larger vegetation cover would be expected naturally, but its dependence on soil type is 
determined by Vs and B functions. The only major exception is soil type 5 for which it 
reduces for 80% and 100% vegetation covers. This can be the result of non-linear stomatal 
control of vegetation to soil moisture, ambient temperature and the PAR (Photosynthetically 
Active Radiation). Such a combined non-linear behavior is not so easy to discern in simple 
theoretical arguments and is beyond the purview of the current focus. The accompanying 
sensible fluxes (not shown) were consistent with the surface energy balance and do not 
have large variations. We naturally expect sensible fluxes to decrease with increasing 

vegetation cover, but this is intertwined into solar energy absorption at the surface through 
the surface albedo, a surf given by: 

“surf = V f“forest +(• - V f )a desert , /m 


while its transpiration characteristics are modulated by the limiting influence of wilting for 



low levels of root zone soil moistures. The interactions lead to a slight increase in sensible 
fluxes with increased vegetation, particularly for soil type 3. It is quite possible that lack of 
the dependence of surface albedo of the soil on the top layer soil moisture is responsible for 
this behavior. 

b) Precipitation: 

The dependence of precipitation on evapotranspiration due to changes in vegetation cover 
and soil type very much affirm the qualitative inferences of Sud et al. (1993 & 1995) and 
Schickedanz and Ackermann (1977). It is shown in the corresponding panels on the right 
hand side of Figure 9a. The simulated precipitation increases with increased vegetation 
cover as well as evapotranspiration. One case of soil type 5 in which the 
evapotranspiration decrease for 80% vegetation cover shows a corresponding decrease in 
precipitation. The spread among the different soil types becomes less as vegetation cover 
increases. As expected, for high vegetation cover, e.g., 99%, both the evaporation and 
precipitation for all soil types tend to be significantly closer. These results suggest that the 
influence of the soil type on evapotranspiration get suppressed as vegetation controls 
influence evaporation. We match of evapotranspiration with observations but the simulated 
precipitation continues to be less than the observed. The error does not seem to be related 
to moisture storage in the atmosphere because the model was shown to simulate realistic 
temperatures and specific humidity profiles (see Section 4.1). Hence, the discrepancy must 
emanate from the prescribed (ARM-CART data) or model bias correcting for the horizontal 
advection (Eq.5, section 2) of moisture. It implicates input data and/or the adjustment 
process itself correcting for SCM biases. Nevertheless, for a given vegetation cover, the 
soils that produce more evapotranspiration also produce more precipitation is abundantly 
evident in these simulations. 

c) Vegetation-Evaporation-Precipitation Relationships 


Here we expand on the main theme of the evapotranspiration-precipitation relationship as ,a 



function of change in vegetation cover and initial soil moisture as well as AP/ AF. First, we 

took the original ensemble of 30 control runs. We plotted the increases in E and P relative 
to 1% vegetation cover for soil type-3 alone as well as the means for all 5-soil types (Figure 

9b). The corresponding AP/AE is also shown in the Figure. The relationships between 
evaporation and rainfall increases as a function of vegetation cover are not so straight 
forward to interpret; nevertheless, AP/AE values remain between 40-60% range, a result 

that is more clearly borne out in the next set of soil moisture controlled simulations. 
Subsequently, we produced the same analysis for 150 simulations in which the initial soil 
moisture was based on the 60% vegetation cover soil moisture initializing runs. The 5- 
ensembles corresponding to five soil moisture initial conditions (described before) were 

again used to produce AP and AE and AP/AE for each of the five ensembles and averaged 

for all the soil types. These provided a very coherent and robust picture (Figure 9c). Since 
the input data, including moisture divergences, are prescribed from the ARM-CART data 
and is identical for each of the 1 50 simulations, the simulations remain heavily constrained 
by the prescribed large-scale flow. Nevertheless, it is modulated by evapotranspiration 
anomalies due to vegetation cover through equations (4a & b) and soil moisture effects. 
One clearly sees how vegetation cover affects the evapotranspiration, precipitation, and the 

ratio AP/AE. Vegetation helps to access the deep soil moisture, but with stomatal controls 

the influence of initial soil moisture on the increase in evapotranspiration are highly non- 
linear. Nevertheless, the relationship of precipitation change in response to 
evapotranspiration change is quite robust. It shows as a positive feedback with about 40- 
60% recovery of evapotranspiration increase through precipitation income. The results are 
quite similar for each soil type (not shown) and even for different initial soil moisture 
values. 

d) Evaporation-Precipitation Relationship 
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Figure 10a shows the 3- week averaged evaporation-precipitation plot for all 180 
simulations. It shows that atmosphere responds to evaporation change in a robust fashion 
regardless of the complexity of the biospheric processes and land-hydrology. This 
response is delightfully simple and linear. It helps us infer that for the Great Plains region, 
roughly 50% of the evapotranspiration is returned back in the form of increased 
precipitation. Indeed regression analysis gave a slope of 0.499, which as a rule of thumb 
can be round off to 50%. We attempted to verify the intercept of 5.06 mm d' 1 precipitation 
dry soil case by running one simulation with zero initial soil moisture. However, the 
simulation produced an evaporation of 0.47 mm d ' through evaporation of rainfall 
interception and soil surface wetting. The corresponding simulated precipitation was 5.26 
mm d 1 (plotted as black dot in Figure 10a). This case too affirms the simple straight line 
relationship. 

e) Moisture distribution in the Vertical 
(i) Vegetation Dependence 

The vertical structure of moisture divergence has some interesting features, Fig. 10b. In 
response to increased vegetation cover, the moisture divergence becomes smaller near the 
surface, i.e., between 960 (surface pressure) and 850 hPa. However, it increases for the 
middle atmosphere, i.e., between 500 and 850 hPa. In other words despite higher 
evapotranspiration with increased vegetation cover, the horizontal moisture divergence 
reduces. The difference is due to enhanced upward transport of moisture by physical 
processes and/or removal by low-level condensation (truly small). An increase in moisture 
divergence in the middle atmosphere represents the influence of moist convection and dry 
turbulent exchanges. These results are naturally complicated by surface albedo decreases 
for higher vegetation cover. Indeed, higher surface energy fluxes would amount to larger 
moist static energy and invigorated moist convection. However, in view of 30 different 
spin ups of initial soil moisture, the land with lower vegetation cover naturally builds up 



higher soil moisture and tries to mitigate evapotranspiration reduction produced by lesser 
vegetation cover. Since, a layer moistened with increased evaporation is expected to 
produce higher moisture divergence, the simulated reduction in moisture divergence 
suggests the dominant role of upward moisture transport by physical processes. In fact not 
only the moist convection would transport the moist air up. but it also brings dry air into 
the near surface layers through convective scale downdrafts. To the extent that the 
moisture loss of lower layer is reflected as moisture gain in the middle layer or even top 
layer (which is miniscule), it is the outcome of the vertical transport. The difference 
between these represents the increase in moisture divergence. Indeed, when we increased 
the depth of the lower layer to 700mb, most of the moisture divergence changes were 
reflected in the lowest layer (not shown). In this way, we note that surface fluxes can 
affect the vertical transport of the moisture through the boundary layer and moist convective 
processes and reveal that despite the evaporation increase, the layer can dry. 

(ii) Evapotranspiration Dependence 

The simulation can be examined also by looking at the dependence of near surface moisture 
divergence on evapotranspiration. Such a result is much easier to interpret. The scatter 
plot in Fig. 11 shows a discernible scatter because the surface energy absorption is a 
function of vegetation cover and that also influences the surface evaporation. As expected, 
the moisture divergence reduces as evapotranspiration increases. This is obviously 
different from its relationship with vegetation cover where it reduced with vegetation cover 
that is responsible for evaporation increase. The upper level moisture divergences, i.e., 
850-200 hPa (not shown), have compensating (opposite sign) effects on the moisture 
divergence simulated in the near surface layers (shown in Fig. 1 1). The sum of the entire 
column moisture divergences yields a smooth graph showing a linear relationship between 
moisture divergence and evapotranspiration. This is also consistent with a linear 
relationship between simulated precipitation and evapotranspiration shown in Fig 10a. 
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5. Discussion and Conclusion. 


The aim of the current investigation is to broadly understand the evaporation-rainfall 
relationship together with the role of soil types, vegetation cover and initial soil moisture. 
If the precipitation were weakly affected by evapotranspiration, the rainfall in the region 
would be largely independent of evapotranspiration. One can intuitively infer that Sahara 
Desert and/or coastal regions with orographic lifting and inland large-scale flow would be 
such regions. But how does the relationship vary in different regions; and how are regions 
with strong evaporation-precipitation feedback distributed around the world and where 
would one have the ability to influence the precipitation by clever management of the 
biosphere remains a puzzle to be addressed and understood. 

All the current single column model simulations for the ARM-CART data confirm a 
positive feedback between evapotranspiration and precipitation. However, there were no 
cases in which precipitation increase exceeded the evapotranspiration increase. Such an 
occurrence would generate a runaway wet environment and evidently early summer season 
of the mid-western United States is not a candidate for such an outcome. We simulate 
larger evapotranspiration in the growing season in response to higher initial vegetation 
cover. The extra water needed for increasing the evapotranspiration would have to be 
supplied through irrigation that draws from external water sources or ground water storage. 
However, because of the evapotranspiration-rainfall relationship and efficient retention of 
rainfall water by the vegetation in the biosphere, the external water supply needs to be 
roughly half of that needed for supplementing evapotranspiration. Obviously, these results 
can not be generalized to other regions. Even for the mid-west July, the new results 
warrant some caution. First we need to establish the robustness of these results with more 
cases. Second, we prescribed the large-scale flow to the SCM. This assumes that the 
regional hydroclimatological changes are not so large as to affect the large-scale. This 
might not hold good on all scales. Nevertheless, we have made an appropriate adjustment 
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for the large-scale forcing of the SCM by introducing a methodology for the grid-scale 
columnar temperature and humidity changes to affect their advective tendencies in the 
horizontal as well as the vertical directions. Consequently, we contend that our study 

addresses a realistic case study in which scale-separation between the large and meso-scales 
has been assumed. 

Indeed, the moisture advected out of the SCM will not a total loss for the United States 
because it has a long way to travel before it leaves the continent. It naturally augments the 
specific humidity of the air mass flowing out of the region and is likely to increase the 
precipitation in the nearby regions (Schickedanz, and Ackermann, 1977). For this reason 
even if the irrigation relies on water that is available a few hundred kilometers around the 
region, the excess rainfall outside the region is likely to recover some of the water usage. 

A natural question is: would the increase in sub-grid scale sensible flux and moist 
processes anomalies due to changes in vegetation cover and evapotranspiration not affect 
the local meso-scale circulation? The answer is: indeed, it would (Pielke and Avissar, 
1990), but the goal of a good parameterization is to handle such effects realistically, 
particularly in the domain-averaged sense. Whether our SCM invoking GEOS II physics 
and McRAS has good enough parameterizations to address such influences remains 
unanswered, but the fact that the SCM responds realistically to several ARM-CART and 
other evaluations provides some credibility to these studies (Sud and Walker, 1999, and 
Ghan et al., 2000). Besides, SSiB and its evapotranspiration sub-model have been 
evaluated thoroughly with several field observations while the new 100-layer land 
hydrology model in invoked to improve the hydrologic accuracy of soil-water conduction. 
Therefore, we are convinced that our evapotranspiration and soil hydrology model is as 
good as any modem state-of-the art SCM can provide. We believe, studies such as this 
help to alleviate some of the confusion created by a variety of answers given by GCMs to a 
typical deforestation such as summarized in Hahmann and Dickinson (1997) and other 
studies of climatic effects of vegetation and/or soil-moisture anomaly scenario. 



In addition to the above limitations of the evaluation procedure and modeling assumptions, 
several inherent parameterization problems continue to hinder an accurate simulation of a 
land-atmosphere interaction. Among the unsolved problems are the influences of subgrid- 
scale effects that include small-scale variability of soil type, soil moisture, as well as the 
natural and agricultural biodiversity, orography and related natural drainage characteristics, 
and surface runoff pathways with its associated effects on the water table and soil moisture 
availability. A highly non-linear response of surface fluxes to such variations makes the 
problem of parameterizing land processes truly hard-to-address. In addition the intrinsic 
variability of biospheric processes, particularly on climatic time-scales, diminishes the 
positive influences of these advances on the improved climate simulation even though 
interactive biosphere climate simulations have begun to show promise (Zeng et al., 1999), 
Under the circumstances, controlled evaluations such as ours may provide a useful answer. 
We plan to continue to test our new results with additional datasets and regional model 
exercises for the mid-west and other important regions of the world. 
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Figure Legends 


Figure 1. 
Figure 2. 

Figure 3a. 
Figure 3b 

Figure 4. 

Figure 5 
Figure 6. 

Figure 7. 


Schematic representation of convective, stratiform and boundary layer clouds in McRAS 
(Adapted from Sud and Walker, 1999). 

SCM testing of McRAS with ARM-CART datasets for Case 1 and Case 3. Figure 2a througl 
2c show simulation errors for the temperature, specific humidity, and precipitation fo 
the SCM-region. 

Coupling of al 00-layer soil hydrology and SSiB with schematic fluxes. 

Soil types in ISLSCP Initiative one data depicted in the three-component mixture viewed 
as soil triangle. Typical mixing ratios of the soils are located in the triangle. 

Influence of vegetation cover on the hydrology and surface fluxes simulated for 1987 
with the ISLSCP Initiative 1 data of the Oklahoma region (a surrogate for ARM-CART 
site). Fig 4a shows actual fields and Fig 4b shows vegetation anomaly minus Control 
differences. 

Same as Fig. 4 except for Sahel region. 

Thirty Initial-time soil-moisture profiles for all case studies. These soil-moistures are 
produced by running each of the five soil types with each of the six vegetation covers for 
six months: January 1 through July 1, with ISLSCP Initiative I forcing data for 1987. 

Simulated error time-series for the Control case: Soil Type 3 and vegetation cover of 
60% with corrections for the model errors: (a) Temperature fields; (b) humidity 
fields; and (c) simulated vis a vis observed precipitation. 
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Figure 8. 


Figure 9a. 


Figure 9b. 


Figure 9c. 


Figure 10a. 


Figure 10b 


Figures 11 


soil Type 3 simulated vis a vis observed fluxes ol (a) evapotranspi, alien and (b) 
sensible heal. The corresponding precipiiation is in 8c. The three lines are (or 1%, 

60%, and 99% vegetation covers versus observailons, which correspond to roughly 
60% vegetation cover. 

Simulated evaporation and precipiiation for the 3-week period for each of the six 
vegetation types using (top) different soil moisture spin up profiles, and (bottom) 60% 
vegetation cover Initial Soil Moisture Profile. 

Same as figure 9b except for a comparison of Soil Type 3 with that of the average of all 
soil types. In this figure the initial soil moisture profiles were the ones from Fig. 6. 

Evaporation and precipitation increases as a function of vegetation fraction with 
prescribed Initial soil moisture profiles. Averaged results for all five soil types are 
used to make the plots. Ap/AE is plotted on top and shows a very robust 40-6 0% 
recovery of increased evaporation through precipitation. 

Evaporation Precipitation Relationship for an ensemble of 5-sets in which each set has 
5 soil types, six vegetation covers and an initial soil moisture profile from 60% 
vegetation-cover initialization produced with a 6-month spinup shown in Fig. 6. 

Averaged moisture divergence between (i) surface and 850 hPa, (ii)850 and 500 hPa, 

(in) 500-200 hPa, and (iv) entire column values as a function of vegetation cover for 
simulation with initial soil moisture profiles of Fig. 6. 

Scatter plots of (i) surface-850 hPa level moisture divergence and (ii) total column 

moisture divergence as a function of surfaoe evaporation regardless of the vegetation 
fraction or initial soil moisture. 



Figure 1 : Schematic Representation of Different Clouds in the GCM 
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Figure 2c. SCM Simulated and Observed Precipitation (mm day"’) for ^ case studies. 
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Evaporation Precipitation 
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Simulated Increase in Moisture Fields 





Simulated Increase in Moisture Fields 



Precipitation (mm/day) 
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Simulated Increase in Moisture Divergence 
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Moisture Divergence (mm/day) 




